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ABSTRACT 

Context. For the past 3 years we have been monitoring 62 K giant stars using precise stellar radial velocity (RV) 
measurements with the 2m Alfred Jensch Telescope of the Thiiringer Landessternwarte Tautenburg (TLS) . 
Aims. To probe the dependence of planet formation on stellar mass by finding planets around intermediate-mass giant 
p. ' stars. 

Methods. We present high accuracy RV measurements of the K1.5 III star 42 Dra and the K4 III star HD 139357. The 
wavelength reference for the RV measurements was provided by an iodine absorption cell placed in the optical path of 
(-H ' the spectrograph. 

Results. Our measurements reveal that the time series of the radial velocity of 42 Dra shows a periodic variation of 
I , 479.1 days with a semiamplitude of K = 112.5 ms" 1 . An orbital solution yields a mass function of f(m) = (5.29 ± 

0.62) x 10 -s solar masses (Mq) and an eccentricity of e = 0.38 ± 0.06. From our template spectra, taken without the 
iodine cell, we determine a metallicity of —0.46 ± 0.04 dex and a stellar mass of 0.98 ± 0.06 M© for this star. 
HD 139357 shows periodic RV variations of 1125.7 days with a semiamplitude K — 159.9 ms -1 . An orbital solution 
yields an eccentricity, e = 0.10 ± 0.02 and mass function, f(m) = (4.79 ± 0.57) x 10" 7 M Q . An iron abundance of 
—0.13 ± 0.04 dex is obtained, and a stellar mass of 1.31 ± 0.24 M© for the parent star is derived. An analysis of the 
HIPPARCOS photometry as well as our Ha core flux measurements reveal no variability with the radial velocity 
period. Keplerian motion is the most likely explanation for the observed radial velocity variations for these stars. 
Conclusions. The K giant stars 42 Dra and HD 139357 host extrasolar planets with "minimum masses" of 3.88 ± 0.85 

■ Jupiter masses Mj up and 9.76 ± 2.15 Mj up , respectively. 

Key words, star: general - stars: variable - stars: individual: 42 Dra, HD 139357 - techniques: radial velocities - stars: 
CO ' late-type - planetary systems 
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0^ . 1- Introduction interpreted as an evidence for MS pollution (Pasquini et al. 

O ■ 2007). 

•• • Until now more than 250 extrasolar planets around solar- Hatzes & Cochran (1993) found first indications of sub- 

I type mam-sequence (MS) host stars have been detected via stellar com p anions ar0U nd giants. They discovered long- 

^ . the RV method. Out of these, the number of discovered period RV var i at ions in three K giants and they proposed 

5^ ■ planetary companions around giant stars is still limited. twQ viable hypotheses for these variations in radial velocity: 

CG ; Our motivation was to enlarge the regime of stellar masses su b s tellar companions or rotational modulation. Recently 

and characteristics of planet host stars by hunting for plan- Hatzeg et al ( 2 ()06) and Reffert at al. (2006) confirmed 

ets around K giants to provide us with important clues to that the imtial RV var i at i ns found by Hatzes & Cochran 

the process of planet formation in a still poorly investigated (1993) in p Gem were in fact due to a p i arie tary companion, 

star domain. This is possible because giant host stars tend Since the first uneq uivocal discovery of the first extrasolar 

to be more massive than the solar-mass objects that have planet around the K giant HD 137759 (1 Dra) by Frink et 

been the traditional targets of most planet searches, thus aL ( 2 002), several giant stars have been found to host gi- 

we can probe the dependence of stellar mass on planet for- ant planetg ( Se tiawan et al. 2003a, 2003b, 2005; Sato al. 

mation. In addition these stars are m a different evolution- 2 003; Hatzes et al. 2005). Currently, a number of groups 

ary status and have different radii and internal structures are act i ve l y searching for extrasolar planets around giant 

compared to their MS counterparts. Finally, in contrast to stars and this has resu i ted in a burst f recent discoveries 

MS stars, the chemical composition of planet host giants ( Sato et aL ( 2 007); Dollinger et al. (2007); Niedzielski et al. 

shows no preference for metal-rich systems, which has been ( 20 Q7); Johnson et al. (2007b); Sato et al. (2008). 



Send offprint requests to: Michaela P. Dollinger e-mail: 

mdoeiiin@eso.org.de 2. Observations and data analysis 

* Based on observations obtained at the 2m Alfred Jensch tele- 
scope at the Thiiringer Landessternwarte Tautenburg 42 Dra and HD 139357 belong to a star sample observed 

since February 2004 from the Thiiringer Landessternwarte 
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Tautenburg (TLS) as part of the Tautenburg Observatory 
Planet Search Programme (TOPS). This programme 
uses the coude echelle spectrograph mounted on the 
2m Alfred Jcnsch telescope. This instrument provides a 
resolving power of R = 67,000 and a wavelength coverage 
of 4700-7400 Angstroms (As) using the so-called "visual" 
(VIS) grism mode. The exposure time for both stars 
ranged between 5-10 minutes depending on the weather 
conditions and this resulted in a signal-to-noise (S/N) ratio 
typically greater than 150. Standard CCD data reduction 
(bias-subtraction, flat-fielding and spectral extraction) 
was performed using standard IRAF routines. An iodine 
absorption cell placed in the optical path provided the 
wavelength reference for the velocity measurements. 

The RVs were computed using a programme that 
largely follows the prescription of Butler et al. (1996). It 
models an observation of the star observed through the cell 
using a high resolution template spectrum of molecular 
iodine taken with a Fourier Transform Spectrometer 
(R « 10 6 ) and a template spectrum of the star taken 
without the iodine cell. The programme takes into account 
possible changes of the instrumental profile (IP) of the 
spectrograph by using the procedure outlined in Valenti et 
al. (1995). Since the IP can also change spatially along a 
spectral order, these were divided into segments (so-called 
spectral chunks). About 130 chunks were used in the final 
analysis. Thus the spatial (and temporal) variations of the 
IP, which can introduce significant RV errors, were mod- 
eled independently for each chunk. Relative wavelength 
shifts between the iodine and stellar template spectra were 
computed using a version of the Fahlman & Glaspey (1973) 
shift-detection algorithm. The RV measurements from all 
chunks were then combined weighted by the inverse square 
of the RV standard deviation for each chunk. The internal 
velocity error of a spectrum is thus the dispersion of all 
segments used for the analysis. We should note that the 
measured Doppler shifts are relative to the stellar template 
and are thus not absolute radial velocities. For the bright 
K giant stars in our programme we typically achieve a 
RV precision of about 3-5 ms _1 . This value refers to the 
internal velocity error. 

To exclude rotational modulation as the cause of the 
RV variations in the 42 Dra and HD 139357 data, we in- 
vestigated stellar activity in these stars via the Ha spectral 
line. Several studies have shown that the central core of Ha 
forms in the chromosphere. Consequently stars of different 
chromospheric activity show a different shape (depth and 
breadth) of the Ha core (Pasquini & Pallavicini 1991). It 
has been widely demonstrated that the Ha line can be used 
as a good "activity indicator" of chromospheric emission 
(Herbig 1985; Pasquini & Pallavicini 1991; Freire Ferrero 
et al. 2004). By measuring the variations of the core of Ha 
with respect to the continuum it is therefore possible to 
investigate the presence and variability of chromospheric 
active regions. 

To determine the chromospheric contribution to the Ha 
profile one has to exclude possible contamination from tel- 
luric lines. The normal procedure is to observe so-called 
"telluric standards" in order to divide out the telluric lines. 
These are typically early-type, rapidly rotating stars which 
are observed every night. However, with the spectra of the 
Tautenburg star sample this subtraction was not done in 



this way because early-type stars have very strong balmcr 
lines so using them for taking out telluric lines is not so 
appropriate. We therefore restricted our Ha measurements 
to the region within ± 0.6 A of the line center. Two ad- 
ditional spectral regions located at ± 50 A provided the 
"continuum" measurement. All of these regions were free of 
telluric contamination. The Ha core measurement was done 
using a programme that shifts all the spectra to a common 
rest wavelength and measures the flux ratio of the defined 
ar c as (Biazzo 2007, private communication) . After running 
the programme, 74 Dra shows a variability of around 1 % in 
Ha activity. Taking into account that this star was not ob- 
served every night in all observing runs this value can only 
give a hint to the internal error of the measurements. Due 
to the Ha activity of only around 1 % and RV variations 
at a very low level, which corresponds to the RV precision 
limit of about 3-5 ms" 1 , 74 Dra was used as an internal 
Ha and RV standard. The RV plot of this star is shown in 
Fig. 1. 

The time analysis of the measured Ha activity will be 
discussed in more detail for both TLS planet candidates in 
§3.1 and §3.2. 



3. Properties of the stars 42 Dra and HD 139357 
and of their planetary companions 

3.1. 42 Dra 

42 Dra (= HD 170693 = HR 6945 = HIP 90344) has a visual 
magnitude of V = 4.83 mag and is classified in SIM BAD 
as K1.5 III star. As for the HIPPARCOS parallaxes, we 
adopt the values derived by van Leeuwen (2007). The cor- 
responding value is 10.36 ± 0.20 mas and this implies an 
absolute magnitude My = —0.09 ± 0.04 mag. As demon- 
strated by van Leeuwen & Fantino (2005) and van Leeuwen 
(2007), this new reduction has substantially improved the 
results for the parallaxes of bright stars, with respect to the 
original values in the ESA catalogue (1997). In the specific 
case of 42 Dra and HD 139357, although the mean paral- 
lax values are close to their previous values, their standard 
errors are now reduced by factors of about 2. 
The stellar parameters of 42 Dra, summarized in Tab. 1 and 
shown in Fig. 2, were cither obtained from the literature or 
derived from our analysis of the stellar spectra taken with- 
out the iodine cell. These high-quality templates allowed 
us to determine accurate Fe abundances [Fe/H], effective 
temperatures T e ff, logarithmic surface gravities logg, and 
microturbulcncc velocities £. 

The results of this analysis for 42 Dra, HD 139357, and 
the rest of the Tautenburg sample as well as a comparison 
with previous studies will be presented in more detail in a 
forthcoming paper. 

Metallicity, T c ff and absolute V-band magnitude, as 
derived from HIPPARCOS parallaxes, were used as 
input values to estimate the mass, age, and radius of each 
programme star by comparing these values to theoreti- 
cal isochrones using a modified version of J^rgensen & 
Lindcgrcn's (2005) method. The result of these method 
is the total Probability Distribution Function (PDF) for 
each stellar property (see Fig. 2). The peak of the PDF 
represents the most probable value of a stellar parameter. 
A detailed description of the procedure is given in da Silva 
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Table 1. Stellar parameters of 42 Dra. The properties of 
the host star are listed in detail. 



Table 2. Radial velocity measurements for 42 Dra. 



Spectral type 


K1.5III 


HIPPARCOS 


mv 


4.833 ± 0.005 


[mag] 




n no -i- n o/i 


[mag] 


B - V 


1.19 ± 0.005 


[maEfl 


Parallax 


10.36 ± 0.20 


[mas] 


Distance 


96.5 ±1.9 


[PC] 


M M 


0.98 ±0.05 


[M ] 


Ri a) 


22.03 ± 1.00 


[Ro] 


Age (a) 


9.49 ± 1.76 


[Gyr] 


rp(a) 


4200 ± 70 


[K] 


[Fe/H] (a) 


-0.46 ±0.05 


[dex] 


logff (a) 


1.71 ±0.05 


[dex] 


micro turbulence' 11 ^ 


1.6 ±0.8 


[kms^ 1 ] 



(a) Dollinger (2008), Dollinger (2009), in preparation 



et al. (2006) Q 

The time series of our RV measurements (see Tab. 2) 
of 42 Dra including the orbital fit is plotted in the upper 
part of Fig. 3. The RV curve of the K giant 42 Dra shows 
an obvious sinusoidal variation. With the exception of one 
point (lower right), possibly caused by stellar oscillations, 
the orbital solution fits the data well. The Lomb-Scargle 
periodogram of the RV measurements is shown in Fig. 4 
and this confirms the presence of a statistically significant 
peak at a frequency of v — 0.00196 cd _1 (P = 510.2 days). 

The "False Alarm Probability" (FAP) of the 510-day 
period peak using the prescription in Scargle (1982) is esti- 
mated to be « 2 x 10" 7 . The FAP was also estimated using 
a "bootstrap randomization technique" (see Kiirster et al. 
1999). The RV values were randomly shuffled keeping the 
times fixed and a Lomb-Scargle periodogram calculated for 
each random data set over the same frequency range shown 
in Fig. 4. After 200,000 "shuffles" there was no instance of 
the random periodogram having power larger than the real 
data set. This confirms that the periodic signal is not due 
to noise or data sampling. 

An orbital solution to the data was made using the 
general non-linear least squares programme GaussFit 
(McArthur et al. 1994). This resulted in a shorter period 
of 479.1 ±6.2 days. The discrepancy between the periods 
is due to the fact that the Lomb-Scargle programme uses 
sine functions (i.e. a circular orbit) to find periodic sig- 
nals in the data. The orbital eccentricity of 42 Dra b is 
e = 0.38 ± 0.06, so a pure sine wave is not the best fit 
to the data. Consequently, the Lomb-Scargle periodogram 
should yield a less reliable period. 

All orbital elements are listed in Tab. 3. The correspond- 
ing mass function is /(m) = (5.29 ± 0.62) x 10~ 8 M Q . 
Using our derived stellar mass of 0.98 ± 0.05 M we calcu- 
lated a "minimum mass" of m sin i = 3.88 ± 0.85 Mj up for 
the companion. The line in Fig. 3 shows the orbital solution 
to the RVs. Fig. 5 shows the phase-folded RV variations and 
orbital fit. 

The RV residuals of 42 Dra are plotted in the lower 
panel of Fig. 3. A periodogram analysis of the residual 
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.672927 


-193 
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3. 
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.380979 


49. 


.9884 


3. 


23 


2454313. 


.374845 


86. 


.8409 


3. 


47 


2454330. 


.329069 


101. 


.4740 


3. 


60 


2454337. 


.334089 


86. 


.0530 


3. 


.15 



Table 3. Summary of all orbital parameters for the plane- 
tary companion to 42 Dra. 



Period [days] 

Tperiastron [JD] 

^[ms^ 1 ] 
a(0-C))[ms- 

e 

w[deg] 

/M[m q ] 

q[AU] 



479.1 ±6.2 
2452757.4 ± 3.7 
110.5 ±7.0 
26.0 

0.38 ±0.06 
218.7 ± 10.6 
(5.29 ±0.62) x 10~ 8 
1.19 ±0.01 



The present implementation of the PDF method is publicly 
available via the web interface 
http : //stev. oapd. inaf . it/param 



RV data showed no additional significant frequencies in 
the data. To check, if rotational modulation might be 
the cause of the observed RV variations we analyzed the 
HIPPARCOS photometry and the Ha activity of 42 Dra. 
Figs. 6 and 7 show the periodogram of the photometry and 
the Ha activity, respectively. Neither the HIPPARCOS 
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photometry nor the Ha data show significant variations at 
the orbital frequency. From the projected rotational veloc- 
ity v sin i < 1.0 kms -1 (de Medeiros et al. 1996) and 
the adopted stellar radius, listed in Tab. 1, we have esti- 
mated the lower limit of the rotational period. This calcu- 
lated value of 1113 days is completely different from the 
orbital period of the planetary companion (see Tab. 3) and 
confirms that the RV variations are not due to rotational 
modulation. The typical lifetime of active zones for K gi- 
ants in general can range between hours to around several 
hundreds of days. 

We note that the sigma of the RVs around the orbital 
solution is about 26 ms" 1 , significantly higher than our RV 
precision. There is also one point which differs significantly 
from the predicted RV of the orbital solution (lower right 
in Fig. 3, near the third RV minimum). This scatter most 
likely arises from stellar oscillations. It is well known that K 
giant stars show RV variations due to p-mode oscillations 
with periods ranging from hours to days (e.g. Hatzes & 
Cochran 1995; de Ridder et al. 2006; Hatzes & Zechmcistcr 
2007). We can use the scaling relations of Kjeldsen & 
Bedding (1995) to estimate the expected amplitude of such 
variations. The velocity amplitude for stellar oscillations 
is predicted to be v osc = ((L/L )/(M/M ))x23.4 cms" 1 
according to Kjeldsen & Bedding (1995) Eq. 7. Using the 
luminosity, calculated with the spectroscopic effective tem- 
perature T e fi, the stellar radius R», and the stellar mass M» 
for 42 Dra results in an amplitude of 22.9 m s _1 , consistent 
with our rms scatter. The characteristic periods for such 
oscillations are expected to be about 1.7 days according to 
the Kjeldsen & Bedding scaling relations. Thus, the scatter 
of the RV data about the orbital solution can be accounted 
for by stellar oscillations. 

3.2. HD 139357 

Another exoplanct orbits the slightly metal-poor giant 
HD 139357 (= HR 5811 = HIP 76311) having the longest 
period of all Tautenburg planet candidates within this K 
giant star sample with a period comparable with the du- 
ration of our survey The stellar parameters of this star are 
summarized in Tab. 4. The PDF diagrams of the stellar 
parameters are shown in Fig. 8. 

The time scries of the RV measurements for HD 139357 
is shown in the upper panel of Fig. 9. The lower panel shows 
the RV residuals after subtracting the orbital solution (see 
below) . 

The radial velocity data reveal a nearly sine-like long- 
period variation. A Lomb-Scargle periodogram of the RV 
data (see Tab. 5) shows significant power corresponding to 
a period of 1124 days. The FAP of this peak is « 10~ 9 . We 
do not show the Lomb-Scargle periodogram as it is readily 
apparent from Fig. 9 that there are periodic RV variations 
in this star. 

An orbital solution yielded a revised period of 1125.7 
days and a slightly eccentric orbit, e = 0.10 ± 0.02. All the 
orbital elements are listed in Tab. 6 and the orbital solution 
is shown as a line in Fig. 9. The mass function (see Tab. 6) 
delivers a "minimum mass" of m sin i = 9.76 ± 2.15 Mj up 
for the substellar companion. It appears that our data win- 
dow covered slightly more than one orbital period. 

The orbital fit to the data is again quite good with an 
rms scatter of only 14 ms _1 . The Kjeldsen & Bedding re- 
lations predict an RV amplitude for stellar oscillations of 



Table 4. The values of the stellar properties of the planet 
host star HD 139357 are listed. 



Spectral type 


K4III 


HIPPARCOS 


my 


5.977 ±0.005 


[mag] 


My 


0.61 ± 0.08 


[mag] 


R T/ 


l.z it u.UUo 


[mag] 


Parallax 


O.40 ± U.oU 


[mas] 


Distance 


118.2 ±4.4 


[pel 

U J 


Mi a) 


1.35 ± 0.24 


[M ] 


Ri a) 


11.47 ±0.75 


[R©] 


Age (a) 


3.07 ± 1.47 


[Gyr] 


rp(a) 


4700 ± 70 


[K] 


[Fe/H] (a) 


-0.13 ±0.05 


[dex] 


lo-gffW 


2.9 ±0.15 


[dex] 


micro turbulence ( - a - > 


1.6 ±0.8 


[kms- 1 ] 



(a) Dollinger (2008), Dollinger (2009), in preparation 

13.7 ms _1 for this star, consistent with the rms scatter 
about the orbit. The Lomb-Scargle periodogram after re- 
moval of the orbital contribution due to the exoplanet shows 
no significant periodic signals. 

To test if the RV period could arise from rotational 
modulation we searched for significant frequencies in the 
106 data points of the HIPPARCOS photometry. Fig. 10 
shows the Lomb-Scargle periodogram of the photometry. 
There is no significant power at the RV frequency. 

Fig. 11 shows the Lomb-Scargle periodogram of the Ha 
variations. Again, there is no significant power near the RV 
frequency. 

4. Discussion 

Our radial velocity measurements indicate that the K 
giant stars 42 Dra and HD 139357 host extrasolar giant 
planets. In the case of 42 Dra, the planet orbit is mod- 
estly eccentric which may argue in favour of Kcplerian 
motion. However, it is not highly eccentric as it is the 
case for i Dra b. It is still possible for surface structure to 
produce such a RV curve. The eccentricity argument for 
HD 139357 is less effective due to its nearly circular orbit. 
Nevertheless, the lack of variability in the HIPPARCOS 
and the Ha data for both stars is more consistent with 
the planet hypothesis. We caution the reader that the 
HIPPARCOS photometry was not simultaneous with 
our RV measurements. Thus, we cannot exclude that 
spots were not present when HIPPARCOS observed 
these stars, but are present now and are causing the RV 
variations. Anyhow, our Ha measurements were made 
simultaneously to the RV data. We therefore believe that 
the detected periods of several hundred days in both 
stars are not due to rotational modulation, but rather to 
planetary companions. 

In the search for exoplanets it is important to investigate, 
if there are differences between discovered planet-hosting 
K giants and MS stars and how this different behaviour of 
MS and giant stars is correlated with the process of planet 
formation. 

The first main difference between MS stars and giants is, 
that giants on average have a somewhat higher mass than 
MS stars surveyed for exoplanets. Taking into account 
that the mass determination in giant stars suffers large 
uncertainties, the masses of planet-hosting giants vary 
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Table 5. Radial velocity measurements for HD 139357. 
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between ~ 0.9 and ~ 3 M Q (da Silva et al. 2006; Dollinger 
2008), while those of MS stars are between ~ 0.75 and ~ 
1.5 M©. 

Another difference between MS and giant stars is their 
large radius in comparison to MS stars. To study this point 
we used the stars of the TLS programme and the previous 
southern giant survey carried out with the Fiber-fed 
Extended Range Optical Spectrograph (FEROS) by 
Setiawan et al. (2004a, 2005). The FEROS and TLS 
giants have radii on average about 10 times larger than 
solar stars (da Silva et al. 2006; Dollinger et al. 2008b). 
High metallicity could favour the migration of the planets 
towards short-period systems; in this case, metal-rich 
stars would have many more short-period planets than 
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Table 6. All orbital parameters for the companion to the 
K giant host star HD 139357. 

Period[days] 1125.7 ±9.0 == 

T perias t r on[JD] 2452466.7 ± 3.2 

Kims- 1 ] 161.2 ±3.2 

a(0 - C))[ms- 1 ] 14.14 

e 0.10 ±0.02 

Lj[deg] 235.4 ±10.6 

/(m)[M Q ] (4.79 ± 0.57) x 10~ 7 

q[AU] 2.36 ±0.02 



metal-poor stars. These planets would be detected among 
the MS stars, but no longer around the giants because 
the star, enlarging its radius as it evolves, would have 
swallowed up them. It was noted in Pasquini et al. (2007) 
(see their Fig. 1) that the metallicity distribution is very 
different for giants than for dwarfs hosting planets with 
comparable long orbital periods, that is, excluding those 
with short orbits. This would indicate that the effects of 
migration, even if present, cannot explain the observed 
trend. 

The third difference has been derived from our high reso- 
lution spectra. The abundance analysis for 42 Dra shows 
that this star is the most metal-poor of our Tautenburg 
sample with [Fc/H] = —0.46 ± 0.05 dex. However, it is 
less metal-poor than HD 13189 (Hatzes et al. 2005) with 
a value of —0.58 ± 0.05 dex and also less metal-poor 
than HD 47536 (Setiawan et al. 2003a) with a value of 
—0.68 ± 0.05 dex. This is interesting because MS stars 
hosting exoplanets tend to be metal-rich compared to 
stars that do not posses exoplanets (Santos et al. 2004). 
Other authors (Schulcr et al. 2005; da Silva et al. 2006) 
have also found evidence that planet-hosting giant stars 
are metal-poor. Pasquini et al. (2007) have shown that 
this characteristic holds for most giants hosting planets. 
We also note that Cochran et al. (2007) found a giant 
planet around the MS star HD 155358 with a metallicity of 
[Fe/H] = —0.68 dex. The basic question is if the increased 
metallicity of planet hosting MS stars enhances planet 
formation, or whether the high metallicity is caused by the 
presence of a planetary system. The correlation with the 
metallicity can thus be either cause or effect. The observed 
metallicity distribution of both types of host stars of 
extrasolar planetary systems should be thus very helpful 
to choose between the two most popular planet formation 
mechanisms: core accretion and gravitational instability. 
In the first case, favoured by the core accretion scenario 
(Pollack et al. 1996), the stars should be overmetallic 
down to their center. Results on MS stars obtained by 
several groups (Fischer & Valenti 2005; Ecuvillon ct al. 
2006) favour this primordial scenario, where stars are 
born in high metallicity clouds. In the second case the 
higher values for the metallicity are not primordial, but 
possibly due to the pollution of the convective envelope 
of the star by the infall of planets and/or planetcsimals 
as well as comets or asteroids (Gonzalez 1998; Laughlin & 
Adams 1997; Gonzalez 2001; Murray et al. 1998; Quillen 
& Holman 2000) onto the star. This pollution can be also 
caused by the total inward migration of a planet onto the 
parent star as well as the transfer of material from the 
disc to the host star as a result of the migration process 
(Goldreich & Tremaine 1980; Lin et al. 1996) or to the 
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break-up and infall of one or more planets in multiple 
systems onto the surface of the star due to gravitational 
interactions (Rasio & Ford 1996) with other companions. 
If the star was polluted by the debris of the planetary 
system, then only the external layers of the atmosphere of 
the star were affected by this pollution (Laughlin & Adams 
1997). Assuming this scenario, the difference in the mass 
of the convective envelope between MS stars and giants 
would explain why the metal excess observed in MS stars 
is not observed among evolved stars. The reason is that 
the metal excess produced by this pollution, while visible 
in the thin atmospheres of solar- type stars, is completely 
diluted in the extended, massive atmospheres of the giants. 
The hypothesis of pollution in combination with the 
convective envelope of evolved stars is a possible explana- 
tion for the different metallicities between MS and giant 
stars, but other reasons like a scenario which includes 
strong differences in planet formation with stellar mass 
and possibly planet migration is plausible and cannot be 
excluded. 



5. Conclusions 

We have found hints that the K giant stars 42 Dra and HD 
139357 host extrasolar planets. Both host stars are metal- 
poor and seem to indicate, in contrast to what is observed 
among MS stars, that giant planets around giant stars do 
not favour metal-rich stars. 
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2 IDL (Interactive Data Language) is a data visualization and 
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Fig. 1. RV measurements of the target star 74 Dra. Given a strict RV limit of < 10 ms 1 this star is per definition 
"constant" and used as RV standard. 
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Fig. 2. Diagram of the total Probability Distribution Function (PDF) of the stellar parameters for 42 Dra. The peak of 
the PDF represents the most probable value of each stellar property. 
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Fig. 3. Radial velocity measurements for 42 Dra (top). The solid line is the orbital solution. The RV residuals (lower) 
seem to contain short-period variations. However the Lomb-Scargle periodogram shows no significant frequency. 
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Fig. 6. Lomb-Scargle periodogram of the HIPPARCOS photometry for 42 Dra. 
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Fig. 7. Lomb-Scargle periodogram of the 42 Dra Ha variations. 
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HD139357 [Fe/H] = -0.13 
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Fig. 8. Stellar parameters of HD 139357. The properties of the host star are listed in detail. 
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Fig. 9. Radial velocity measurements and RV residuals for HD 139357. 
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Fig. 10. Lomb-Scargle periodogram of the HIPPARCOS photometry for HD 139357. 
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Fig. 11. Lomb-Scargle periodogram of the Ha variations for HD 139357. 



